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摘要 : ATF) UMRAO 数据 库 的 22 个 焰 变 体 源 的 偏振 数据 ， 根 据 Lazarian & Pogosyan 
的 理论 模型 进行 分 析 研 究 了 炮 变 体 射 电波 段 的 偏振 随 波 长 的 变化 ,得 出 不 同类 的 灯 变 体 偏振 
随 波 长 的 变化 的 原因 。 其 结果 表明 ， 与 Lazarian & Pogosyan 的 理论 模型 吻合 较 好 。 由 此 我 
们 得 出 如 下 结论 :〈1) 反常 去 偏振 焰 变 体 源 中 反常 法 拉 第 旋光 起 主导 作用 ,反常 去 偏振 相对 
较 少 ; (2) 热 辐射 与 同步 辐射 混合 作用 ， 高 频 波段 的 偏振 度 就 小 于 低频 波段 的 偏振 度 《 有 反常 
去 偏振 );〈(3) 当 被 观测 波段 波长 小 于 吸 积 盘 热 辐射 所 对 应 的 维 因 波长 时 ， 高 频 波段 的 偏振 
度 就 大 于 低频 波段 的 偏振 度 〈 和 常规 去 偏振 )。 
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作为 活动 星系 核 HEAR) 的 一 个 子 类 ,辐射 与 观察 者 的 视线 相 一 致 的 次 变 体 CBlazars ) 
被 相对 论 喷 流 所 主导 叫 。 具 有 同步 辐射 的 相对 论 顺 流 和 强烈 的 线 偏振 。 线 偏振 度 可 以 计算 为 


p = -全 ， 其 中 a 是 同步 辐射 的 光谱 指数 ， 磁 场 分 布 可 以 稍微 的 改变 偏振 中。 在 本 文中 ， 利 
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用 UMRAO (https://dept.astro.lsa.umich.edu/datasets/umrao.php〉 数据库 给 出 的 观测 结果 ， 研 
F blazar 的 偏振 度 与 波长 的 变化 。 UMRAO 数据 库 提 供 了 一 个 很 好 的 研究 射 blazar 电波 段 的 
偏振 特性 。 

由 于 法 拉 第 散射 ， 去 偏振 对 同步 辐射 源 的 偏振 度 有 一 定 的 作用 。 偶 振 的 形式 可 以 为 p x 
e-X ， 其 中 4 是 观测 波长 g。 法 拉 第 旋转 量 (RM) 可 能 与 磁场 拓扑 结构 有 关中 中 ， 并 存在 一 
种 特殊 的 磁场 扭曲 效应 [01。 另 外 ， 在 近期 的 研究 中 ，MarscherA.P 等 提出 了 法 拉 第 RM 中 
的 油 流 特性 ， 推 导出 了 一 个 p ce /的 形式 ， 其 中 7 与 汕 流 能 级 联 的 指数 有 关 团 ， 一 些 模拟 
也 表明 ， 满 流 可 能 是 blazar 偏振 变化 的 一 个 原因 外 。 近 年 来 ， 综 合 利 用 统计 相关 函数 研究 了 
三 维 各 同 异 性 磁 流 体 动力 学 MHD)， 来 说 明 满 流 对 偏振 随 波 长 变化 的 影响 ， 也 表明 可 以 应 
用 Lazarian& Pogosyan (2016) 给 出 的 MHD mE ARRIT UMRAO 数据 库 的 射电 波 
段 的 偏振 。 

本 文 是 研究 射电 波段 中 的 去 侦 振 特性 ， 揭 示 偏 振 随 波 长 变化 的 物理 原因 。 我 们 收集 了 
UMRAO 中 2009 年 9 月 份 到 2012 年 5 月份 的 偏振 数据 ,并 通过 px aAX “的 圭 律 拟 合 多 波段 
偏振 度 ， 其 中 b EMES WIERA. RIIE E Lazarian& Pogosyan (2016 F) 
给 出 的 法 拉 第 旋转 的 物理 条 件 ; 然后 ， 给 出 的 与 波长 有 关 的 去 偏振 特征 的 物理 原因 。 


2. 样本 选择 和 拟 合 结果 


一 般 认 为 ， 耀 变 体 的 侦 振 是 由 相对 论 电 子 同 步 辐射 产生 的 ， 有 单 电 子 同 步 辐射 和 集体 电 
子 同 步 辐射 之 分 。 由 于 兆 变 体 的 偏振 是 由 大 量 电子 产生 的 , 所 以 用 集体 电子 同步 辐射 来 研究 
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He Ae ASS) Aiea de AZ Aiea FI o 
HH EE aS HE] SY, BETERIKOA: 
jO) = SV) N(y,0)dydQe (1) 
KF, Og = 2xsin6qd0, N(y,96) 是 电子 能 谱 , SQV) 是 单 相对 论 电 子 同步 辐射 谱 。 通常 取 y 为 y1 < 
y < 加 。 电 子 能 谱 服从 宕 率 分 布 N(y,8) =N, yg), n 为 电子 的 能 谱 指 数 ，g(9) 为 在 9 


方向 上 的 电子 分 布 函数 。 则 集体 电子 同步 辐射 产生 的 偏振 度 为 : 


P= ie (2) 
EMEY) <K v < v(y;)， 则 集体 电子 同步 辐射 的 偏振 度 为 
_ n+l 
n+7/3 


次 变 体 的 连续 谱 通 常 满足 F, x v-*，a 为 谱 指数 ， 所 以 ，n 和 a 的 关系 为 n= 2a+1. 
此 瘤 变 体 的 偏振 度 可 表示 为 口 : 


_ ati 
P= pear (4) 


KENEH, ESE) Ais SE SHS CBR) 无 关 。 
DAE REP ESD, MERS EELS) RONDA. Ballard 等 人 0 提出 了 一 
PERG SEA LY OUR ide, FERIA TUN: 


Í 
P(v) = Hq) POTS (5) 


FOP (vy) 4a Hi APE, FI Le HE RE. C5) 式 中 可 以 看 出 ， 偏 振 度 与 频率 ( 波 
长 》 有 关 ， 偏 振 度 是 随 波 长 变化 而 变化 的 WDP (wavelength dependence polarization )。 

产生 WDP 的 原因 ， 即 产生 Hv) 的 原因 ， 解 释 有 多 种 ， 人 磁场 浦 流 也 会 对 偏振 产生 一 定 的 
影响 (去 偏振 的 内 部 原因 )。 产 生 WDP 行为 的 原因 有 两 种 : CG) 热 辐射 和 同步 辐射 混合 作 
用 ， 使 得 产生 WDP 行为 ; (2) 同步 辐射 中 本 来 就 存在 WDP 行为 0]。 耀 变 体 的 谱 指 数 w 一 
般 在 0 到 2 之 间 ， 对 应 于 〈4) 式 的 偏振 度 为 60% 到 82%。 然 而 实际 测 得 的 射电 波段 的 偏振 
度 一 般 低 于 10%， 则 ICv) 起 到 了 主要 的 作用 。 

Lazarian 和 Pogosyan03 分 别 讨论 了 法 拉 第 旋转 主导 和 消 流 主导 下 对 偶 振 的 影响 , 偏振 随 
波长 的 变化 : (1) 在 法 拉 第 旋转 主导 下 ， 偏 振 随 波长 变化 的 关系 为 DCA) x 和 -下 或 者 pQ) x 
和 0 (2) 在 清流 主导 下 ， 偏 振 随 波长 的 变化 关系 为 pO) x ACT /2 NE pA) x 
和 TO mg 其 中 nt 分别 对 应 于 汕 流 的 能 谱 指数 和 法 拉 第 RM 的 能 谱 指 数 ， 偏 振 随 波 
长 变化 的 关系 服从 窜 率 形式 。 在 文献 [14] 中 用 p= aA 来 拟 合 数据 ， 我 们 仍 采 用 此 形式 来 拟 
A UMRAO 数据 库 中 的 射电 数据 ， 得 到 不 同 源 对 应 的 b 值 ， 再 得 到 与 之 相对 应 的 m 值 。 

如 果 我 们 假定 mg = 0， 则 汝 流 主导 与 法 拉 第 旋转 主导 下 的 偏振 随 波 长 变化 保持 一 
BUPA) x ACI; my > 1 的 效果 与 mg = 1 的 效果 一 致 ， 则 取 myg = 103]。 

在 选取 样本 时 ， 由 于 不 同 探测 器 之 间 的 不 协调 性 ， 我 们 只 选取 了 UMRAO 

(https://dept.astro.lsa.umich.edu/datasets/umrao.php ) 数据 库 的 4.8GHz、8GHz、14.5GHz ‘in 
振 数 据 ， 且 只 选 了 2009 年 9 月 份 到 2012 年 5 月 份 的 偏振 数据 。 在 众多 AGN 源 中 ， 我 们 只 
选取 利于 研究 的 22 个 兆 变 体 源 ， 在 选取 的 22 个 兆 变 体 源 的 数据 中 ， 观 测 的 平均 周期 为 1 
KR, H 2009 年 9 月 份 到 2012 年 5 月 份 没 有 表现 出 激烈 的 爆发 。 由 于 三 个 波段 的 偏振 数据 不 
是 准 同时 性 的 ， 且 为 了 研究 的 方便 ,我们 对 它们 的 偏振 和 误差 分 别 求 平均 值 ， 得 到 3 个 数据 
Rio FERRET RBA, 3 个 数据 点 足以 说 明 暴 率 情 况 ， 拟 合共 有 可 行 性 。UMRAO 数据 库 
的 数据 来 源 于 NOEMA 射电 望远镜 观测 的 数据 ， 其 主要 目的 是 用 于 校准 其 他 射电 望远镜 获 


得 的 数据 ， 且 其 中 的 线 侦 振 为 部 分 线 偏振 ， 数 据 库 仅 包含 每 日 观测 的 侦 振 度 P>3sigmaP 的 
数据 。 


图 1 b<O 时 的 耀 变 体 肾 率 拟 合 情 况 ， 横 轴 为 波长 ， 纵 轴 为 偏振 度 


Fig.l Fitting power factor of blazars source for b<0, where the horizontal axis is the wavelength and the vertical axis is the polarization 


degree. 


由 图 1 可 以 看 出 b 值 小 于 0 偏振 度 随 波 长 而 降低 dp/dà > 0， 表 现 出 反常 去 偏振 。 很 难 
用 MHD 消 流 和 法 拉 第 旋转 简单 的 来 解释 产生 dp/d > 0 的 物理 机 制 。 反 第 的 去 偏振 特性 ， 
即 偏振 随 波 长 变 长 而 增加 ， 在 文献 [14] 的 文章 中 ， 光 学 波段 表现 的 很 少 ， 而 本 文 的 射电 波段 
却 有 相对 较 多 的 表现 ， 大 约 占 炊 变 体 源 的 27%， 接 近 1/3 。 这 种 现象 大 致意 味 着 大 约 有 1/3 
的 次 变 体 喷 流 中 ， 吸 积 盘 的 内 部 的 偏振 比 外 围 的 偏振 低 ， 相 当 于 外 围 磁 场 强 ， 越 徘 内 磁场 越 
弱 ， 这 是 很 反常 的 。 这 种 反常 ， 说 明 很 有 可 能 这 些 1/3 的 光 变 体 的 吸 积 盘 或 噶 流 中 的 磁场 发 
生 了 扭曲 。 这 种 扭曲 使 得 常规 的 同步 辐射 和 MHD 模型 中 的 磁场 在 数据 结果 中 与 预测 的 不 一 
致 ， 完 其 原因 ， 我 们 认为 有 可 能 是 应 动 磁 重 联 造成 的 。 即 导致 了 磁场 拓扑 结构 发 生 改 变 ， 将 
磁 能 转化 为 等 离子 体 的 动能 和 热能 ,并 加 速 等 离子 体 中 的 部 分 带电 粒子 , 吸 积 盘 内 的 磁场 一 
部 分 被 耗 散 掉 ， 耗 散 的 磁场 转化 为 等 离子 体 的 动能 和 热能 ， 并 加 速 等 离子 体 中 的 部 分 粒子 ， 
使 得 盘 内 磁场 减 小 ,小 到 低 于 外 围 磁 场 强 度 。 最 后 形成 热 辐射 与 同步 辐射 相 混 合 ， 进 而 反常 
法 拉 第 旋光 起 主导 作用 。 
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图 2 b>0 时 的 炮 变 体 源 窜 率 拟 合 情 况 ， 模 轴 为 波长 ， 纵 轴 为 偏振 度 


Fig.2 blazars source power law fitting for b>0, where the horizontal axis is the wavelength and the vertical axis is the polarization 


degree. 

图 2 给 出 了 b>0 的 情况 。 有 14 NERA. Hite PE ER MK dp/dà < 0， 表 现 出 常 
规 去 偏振 。 大 致 表明 在 这 些 泡 变 体 喷 流 中 ， 吸 积 盘 的 内 部 的 偏振 比 外 围 的 偏振 高 ， 外 围 磁 场 
弱 ， 越 靠 内 磁场 越 强 ， 服 从 同步 辐射 的 磁场 特性 。 这 种 情况 大 约 占 70%， 属 于 大 多 数 情况 。 
当 观 测 波段 小 于 吸 积 盘 热 辐射 所 对 应 的 维 恩 波 长 〈 热 辐射 的 最 大 波长 ) 时 ,高 频 波 段 的 偏振 
度 就 大 于 低频 波段 的 偏振 度 ， 与 参考 文献 [17] 中 光学 波段 所 预测 的 一 致 。 


201809.00123v1 


chinaXiv 


number 


图 3 b 值 的 分 布 情况 。 


Fig.3 Distribution of b values. 
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图 4 ESR UE b 值 情况 。 其 中 b>0, 代表 dp/ 
dà < 0; b<0, 代 表 dp/da > 0。 


Fig.4 b-values for each blazar source. Where b>0, 


stands for dp/da < 0; b<0, stands for dp/da > 0. 


b value 


source name 


图 3 给 出 了 b 值 的 分 布 直方 图 , 由 图 可 以 看 出 , b>0 的 个 数 多 与 b<0 的 个 数 , 这 表明 dp/ 
dà > 0 在 偏振 随 波 长 的 变化 中 起 了 主导 作用 。b 值 大 致 分 布 在 0 到 0.8 之 间 ， 主 要 集中 分 布 
在 0 到 0.4 之 间 , 图 4 表明 每 个 浴 变 体 源 b 值 得 分 布 情况 ,可 以 明显 看 出 只 有 6 个 源 位 于 b<0， 
其 余 的 都 是 b>0。 


图 6 消 流 磁场 主导 时 消 流 指数 m 的 分 布 情况 ， 每 个 m 值 
对 应 于 相应 的 b 值 ， 注 意 到 m=2/3 对 应 于 三 维 各 向 异性 


Kolmogorov #%. 
Fig.6 Distribution of the MHD turbulent index m in the 
turbulent magnetic field. Each value of m is derived from the 
fitting value of b. We note that m=2/3 corresponds to the 


three-dimensional anisotropic Kolmogorov scaling. 
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图 5 平均 磁场 主导 时 消 流 指数 m 的 分 布 情况 , 每 个 m 值 对 
应 于 相应 的 b 值 ， 注 意 到 m=2/3 对 应 于 三 维 各 向 异性 


表 1 样本 拟 合 结果 与 FDP 特性 总 结 。 


Table 1 Summary of sample fit results and 


Kolmogorov 谱 。 


Fig.5 Distribution of the MHD turbulent index m in the 
. . | FDP characteristics. 

regular magnetic field dominated case. Each value of m is 

derived from the fitting value of b. We note that m=2/3 


corresponds to the three-dimensional anisotropic Kolmogorov 


scaling. 
Name Type Redshift b mR mr 
0133+476 (DA 55) BL Lac 0.859 0.10+0.12 0.79 0.90 
0215+015 (OD 026) BL Lac 1.715 0.37+0.07 0.26 0.63 
0235+164 BL Lac 0.940 -0.47+40.49# 1.94 1.47 
0316+413 (3C 84) Blazar 0.01756 -0.27+40.03* 1.54 1.27 
0528+134 BL Lac 2.070 -0.09+£0.18* 1.18 1.09 
0716+714 BL Lac 0.300 0.38+0.41 0.25 0.62 
0727-115 BL Lac 1.591 0.22+0.29 0.56 0.78 
0851+202 (OJ 287) BL Lac 0.306 0.73+0.35 -0.45 0.27 
1127-145 BL Lac 1.18906 -0.47+0.1 8” 1.94 1.47 
1226+023 (3C 273) BL Lac 0.158 0.34+0.31 0.31 0.66 
1253-055 (3C 279) Blazar 0.53620 -0.40+0.05* 1.79 1.40 
1308+326 BL Lac 0.996 0.03+0.02 0.95 0.97 
1502+106 (OR 103) BL Lac 1.83928 0.53+0.39 -0.06 0.47 
1510-089 Blazar 0.361 -0.16+0.3 1” 1.31 1.16 
1633+382 (4C+38.41) Blazar 1.81313 0.29+0.27 0.41 0.71 
1641+399 (3C 345) Blazar 0.59476 0.05+0.14 0.90 0.95 
1749+096 (OT 081) BL Lac 0.322 0.1940.08 0.62 0.81 
2022-077 BL Lac 1.388000 0.37+0.42 0.25 0.63 
2200+420 (BL LAC) BL Lac 0.069 0.08+0.00 0.84 0.92 
2230+114 (CTA 102) BL Lac 1.037 0.16+0.15 0.68 0.84 


VE: m 表示 在 平均 磁场 占 优 情 况 下 得 到 的 MHD mA, m HAN TET a AC BS BIAS MHD diate Al. CH) 代表 dp/ 
dà > 0 HEAR 


由 图 5 可 看 出 ,平均 磁场 主导 时 消 流 指数 m 的 分 布 情况 ,b 值 所 对 应 的 m 值 分 布 直方 图 。 


由 图 6， 满 流 磁场 主导 时 清流 指数 m 的 分 布 情况 ，b 值 所 对 应 的 m 值 分 布 直方 网 。 看 起 来 
观测 样本 的 拟 合 结果 与 平均 磁场 主导 (法拉第 旋转 占 优 ) 情况 和 消 流 磁场 主导 〈 法 拉 第 旋转 
占 优 ) 情况 大 致 一 致 。 表 1 中 ,我们 列 出 了 每 个 源 偏振 观测 的 结果 特性 ， 发 现 OJ 287 和 OR 
103 的 mx 值 小 于 0， 这 与 文献 [1 和 1 中 的 结果 人 不符， 其 原因 不 尽 相 同 。 可 能 结 末 是 因为 数据 点 
的 不 足 ， 或 者 是 在 数据 的 时 间 段 内 偏振 变化 相对 较 大 ， 从 而 导致 取 数 据 的 平均 值 误差 较 大 。 


3. wit 


在 拟 合 的 结果 中 ， 有 两 个 源 0420-014、0805-077 的 结果 为 p。( 偏 振 随 波 长 变化 不 明显 )。 
有 14 个 源 的 dp/dà < 0, 6 个 源 dp/da > 0。 次 变 体 的 偏振 主要 来 自 同 步 辐射 , 但 星际 介质 、 
吸 积 盘 、 宽 线 区 等 也 会 对 偏振 造成 一 定 的 影响 。 在 本 文中 , 我们 只 考虑 磁场 淇 流 和 法 拉 第 旋 
光 的 影响 ， 并 忽略 了 以 上 效应 。 由 于 只 把 电子 能 谱 作 为 简单 的 容 率 形式 ， 所 以 内 材 偏 振 度 与 
波长 无 关 ， 对 结果 无 影响 。 我们 观察 到 样本 中 dp/dA < 0 的 去 偏振 特性 ， 也 注意 到 表 1 中 列 
出 了 一 些 dp/da > 0 的 情况 , 很 难 应 用 MHD 湛 流 和 法 拉 第 旋转 来 解释 产生 dp /da > 0 的 物 
理 机 制 ， 应 用 反常 的 去 偏振 特性 ， 即 偏振 随 波 长 变 长 而 增加 , 扭曲 的 磁场 可 以 减少 反常 的 法 
拉 第 旋转 效应 3。 因 此 ，dp/dx > 0 的 去 偏振 情况 是 可 以 及 生 的 。 最 近 ， 在 一 些 兆 变 体 喷 流 
中 发 现 了 常规 去 偏振 Cdp/dA<0) 和 反常 去 偏振 Cdp/dA> 0) RENS, ee HE a die HY He 
变 体 源 中 反常 法 拉 第 旋光 起 主导 作用 ,反常 去 偏振 相对 较 少 。 本 文 样品 中 得 到 的 去 偏振 情况 ， 
表明 了 喷 流 中 偏振 的 特性 : 热 辐射 与 同步 辐射 混合 后 , 高 频 波段 的 偏振 度 就 小 于 低频 波段 的 
偏振 度 ( 反 常 去 偏振 )， 当 观测 波段 小 于 吸 积 盘 热 辐射 所 对 应 的 维 恩 波长 ( 热 辐 射 的 最 大 小 
长 )》 时 ， 高 频 波段 的 偏振 度 就 大 于 低频 波段 的 偏振 度 t2 (常规 去 偏振 )。 

Lazarian & Pogosyan 的 理论 分 析 表 明 ， 法 拉 第 旋转 波动 来 自 各 问 异 性 MHD miio R 
们 在 耀 变 体 样品 中 的 拟 合 统计 结 果 表 明 ， 射 电波 段 的 去 偏振 大 致 服从 Kolmogorov 谱 。 去 偏 
振 和 相关 的 潮流 特性 ， 表 明了 不 同 光 变 体 源 的 偏振 多 样 性 。 

对 某 一 耀 变 体 而 言 ， 不 同 波长 的 射电 偏振 并 不 是 在 时 间 上 同步 或 准 同步 观测 得 到 的 ， 那 
么 ， 本 文 的 分 析 和 结果 以 及 物理 依据 不 是 很 合适 的 。 另 外 ， 在 本 文中 可 以 看 到 ,样本 中 大 部 
分 滩 变 体 ， 不 同 波长 的 偏振 度 变 化 幅度 并 不 大 ， 数 据 只 有 三 个 波段 ， 波 段 范 围 也 很 窄 ， 不 利 
于 拟 合 和 分 析 。 

文章 中 提 到 的 有 关 磁 场 潮流 和 磁场 有 序 度 H(v)， 我 们 将 在 接 下 来 的 文章 中 进一步 研究 。 
ICDJ 和 汕 流 对 磁场 起 到 的 具体 作用 及 物理 机 制 ， 有 着 重要 的 研究 价值 。 


TE: 表 一 中 列举 的 源 的 类 型 ， 期 初 由 于 对 SIMBAD 数据 库 的 不 熟悉 ， 从 而 导致 了 对 源 
的 分 类 的 误 判 ， 对 审 稿 人 造成 的 误解 和 麻烦 ， 我们 深 表 蒜 意 。 在 今后 的 学 习 中 ， 我 们 会 熟悉 
数据 库 的 使 用 ， 不 会 再 出 现 如 此 类 似 的 错误 。 由 此 ， 我 们 再 次 深 表 抱歉 。 

致谢 : 感谢 编辑 和 审 稿 人 对 我 们 稿件 “次 变 体 在 射电 波段 的 偏振 研究 ”提供 的 意见 。 这 
些 意见 对 我 们 进行 论文 修改 和 提高 论文 质量 非常 有 帮助 。 我 们 根据 审 稿 人 的 意见 对 论文 做 了 
仔细 修改 。 在 今后 的 学 习 中 ， 我 会 广 读 文献 来 加 强 目 己 知识 点 的 不 足 。 
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PolarizationStudy of Blazar in Radio frequency band 
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Abstract: In this paper, the polarization data of 22 ray sources in UMRAO database are used to 


analyze thechange of polarization with wavelengthin the radiofrequency bandonblazars, according 


to the theoretical model of Lazarian & Pogosyan. The cause of the change ofpolarizationwith 


wavelengthaboutdifferent class of blazars. The results show that it is in good agreement with the 


theoretical model of Lazarian & Pogosyan. From this, we conclude that: (1) the anomalous 


Faraday rotation in the anomalous depolarization source is dominant, and the anomalous 


depolarization is relatively small; (2) when the combination of thermal radiation and synchrotron 
radiation, the degree of polarization in the high frequency band is less than the of 
polarizationdegree in the low frequency band (abnormal depolarization); (3) when the wavelength 
of the observed frequencyband is smaller than the Wien wavelength of the thermal radiation of the 
accretion disk, the polarizationdegree of the high frequency band is greater than the 
polarizationdegree of the low frequency band (conventional depolarization). 
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